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(54) Adaptive equalization and Interpolated timing recovery in a sampled amplitude read channel 
for magnetic recording 

(57) A sampled amplitude read channel for mag- 
netic disk recording which asynchronously samples (24) 
the analog read signal, adaptively equalizes (B103) the 
resulting discrete time sample values according to a tar- 
get partial response, extracts synchronous sample val- 
ues through interpolated timing recovery (8100), and 
detects digital data from the synchronous sample val- 
ues using a Viterbi sequence detector (34) is disclosed. 
With interpolated timing recovery, the equalizer filter 
and its associated latency are removed from the timing 
recovery loop, thereby allowing a higher order discrete 
time filter and a lower order analog filter. 
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Description 

FIELD OF INVENTION 

5 The present invention relates to the control of magnetic disk storage systems for digital computers, particularly to 
a sampled amplitude read channel that employs asynchronous sampling of an analog read signal, adaptive discrete 
time equalization, and interpolated timing recovery. 

CROSS REFERENCE TO RELATED APPLICATIONS AND PATENTS 

10 

This application is related to other co-pending U.S. patent applications, namely application serial numbers 
08/440,515 entitled "Sampled Amplitude Read Channel For Reading User Data and Embedded Servo Data From a 
Magnetic Medium." 08/341.251 entitled "Sanpled Amplitude Read Channel Comprising Sample Estimation Equaliza- 
tion. Defect Scanning, Channel Quality. Digital Servo Demodulation. PID Filter for Timing Recovery, and DC Offset 

15 Control," 08/313,491 entitled "Improved Timing Recovery For Synchronous Partial Response Recording," and 
08/533,797 entitled "Improved Fault Tolerant Sync Mark Detector For Sampled Amplitude Magnetic Recording." This 
application is also related to several U.S. patents, namely U.S. Pat. No. 5.359.631 entitled "Timing Recovery Circuit for 
Synchronous Waveform Sampling." 5,291,499 entitled "Method and /\pparatus for Reduced-Complexity Viterbi-Type 
Sequence Detectors." 5,297,184 entitled "Gain Control Circuit for Synchronous waveform Sampling." 5.329.554 enti- 

20 tied "Digital Pulse Detector." and 5.424.881 entitled "Synchronous Read Channel." All of the above-named patent appli- 
cations and patents are assigned to the same entity, and all are incorporated herein by reference. 

BACKGROUND OF THE INVENTION 

25 In magnetic storage systems for computers, digital data serves to modulate the current in a read/write head coil in 
order to write a sequence of corresponding magnetic flux transitions onto the surface of a magnetic medium in concen- 
tric, radially spaced tracks at a predetermined baud rate. When reading this recorded data, the read/write head again 
passes over the magnetic medium and transduces the magnetic transitions into pulses in an analog read signal that 
alternate in polarity. These pulses are then decoded by read channel circuitry to reproduce the digital data. 

30 Decoding the pulses into a digital sequence can be performed by a simple peak detector in a conventional analog 
read channel or. as in more recent designs, by a discrete time sequence detector in a sampled amplitude read channel. 
Discrete time sequence detectors are preferred over simple analog pulse detectors because they compensate for inter- 
symbol interference (ISI) and are less susceptible to channel noise. As a result, discrete time sequence detectors 
increase the capacity and reliability of the storage system. 

35 There are several well known discrete time sequence detection methods including discrete time pulse detection 
(DPD), partial response (PR) with Viterbi detection, maximum likelihood sequence detection (MLSD), decision-feed- 
back equalization (DFE), enhanced decision-feedback equalization (EDFE). and fixed-delay tree-search with decision- 
feedback (FDTS/DF). 

In conventional peak detection schemes, analog circuitry, responsive to threshold crossing a derivative Infbrma- 

40 tion, detects peaks in the continuous time analog signal generated by the read head. The analog read signal is "seg- 
mented" into bit cell periods and interpreted during these segments of time. The presence of a peak during the bit cell 
period is detected as a "1 " bit. whereas the absence of a peak is detected as a "0" bit. The most common errors in 
detection occur when the bit cells are not correctly aligned with the analog pulse data. Timing recovery, then, adjusts 
the bit cell periods so that the peaks occur in the center of the bit cells on average in order to minimize detection errors. 

45 Since timing information is derived only when peaks are detected, the input data stream Is normally run length limited 
(RLL) to limit the number of consecutive "0" bits. 

As the pulses are packed closer together on the concentric data tracks in tiie effort to increase data density, detec- 
tion errors can also occur due to intersymbol interference, a distortion in the read signal caused by closely spaced over- 
lapping pulses. This interference can cause a peak to shift out of its bit cell, or its magnitude to decrease, resulting in a 

50 detection error. The ISI effect is reduced by decreasing the data density or by employing an encoding scheme that 
ensures a minimum number of "0" bits occur between "1" bits. For example, a (d.k) run length limited (RLL) code con- 
strains to d the minimum number of "0" bits between "1" bits, and to k the maximum number of consecutive "0" bits. A 
typical (1 ,7) RLL 2/3 rate code encodes 8 bit data words into 12 bit codewords to satisfy the (1 ,7) constraint. 

Sampled amplitude detection, such as partial response (PR) with Viterbi detection, allows tor increased data den- 

55 sity by compensating for intersymbol interference and the effect of channel noise. Unlike conventional peak detection 
systems, sampled amplitude recording detects digital data by interpreting, at discrete time Instances, the actual value 
of the pulse data. To this end, the read channel comprises a sampling device for sampling the analog read signal, and 
a timing recovery circuit for synchronizing the samples to the baud rate (code bit rate). Before sampling the pulses, a 
variable gain amplifier adjusts the read signal's amplitude to a nominal value, and a low pass analog filter filters the read 
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signal to attenuate aliasing noise. After sampling, a digital equalizer filter equalizes the sample values according to a 
desired partial response, and a discrete time sequence detector, such as a Viterbi detector, interprets the equalized 
sample values in context to determine a most likely sequence for the digital data (i.e., maximum likelihood sequence 
detection (MLSD)). MLSD takes into account the effect of ISI and channel noise in the detection algorithm, thereby 
5 decreasing the probability of a detection error. This increases the effective signal to noise ratio and. for a given (d.k) 
constraint, allows for significantly higher data density as compared to conventional analog peak detection read chan- 
nels. 

The application of sampled amplitude techniques to digital communication channels is well documented. See Y. 
Kabal and S. Pasupathy, "Partial Response Signaling", IEEE Trans. Commun. Tech., Vol. COM-23. pp.921-934, Sept. 
10 1975; and Edward A. Lee and Davkl G. Messerschmitt. "Digital Communication". Kluwer Academic Publishers, Boston, 
1990; and G.D. Forney, Jr., "The Viterbi Algorithm", Proc. IEEE, Vol. 61, pp. 268-278, March 1973. 

Applying sampled amplitude techniques to magnetic storage systems is also well documented. See Roy 0. Cideci- 
yan. Francois Dolivo. Walter Hirt. and Wolfgang Schott. "A PRML System for Digital Magnetic Recording". IEEE Journal 
on selected Areas in Communications, Vol. 10 No. 1, January 1992. pp.38-56; and Wood et al, "Viterbi Detection of 
15 Class IV Partial Response on a Magnetic Recording Channel IEEE Trans. Commun, Vol. Com-34, No. 5, pp. 454- 
461, May 1986; and Coker Et al, "Implementation of PRML in a Rigid Disk Drive". IEEE Trans, on Magnetics, Vol. 27. 
No. 6, Nov. 1991; and Carley et al, "Adaptive Continous-Time Equalization Followed By FDTS/DF Sequence Detec- 
tion", Digest of The Magnetic Recording Conference, August 15-17, 1994. pp. C3; and Moon et al, "Constrained-Com- 
plexity Equalizer Design for Fixed Delay Tree search with Decision Feedback", IEEE Trans, on Magnetics, Vol. 30, No. 
20 5, Sept. 1994; and Abbott et al. "Timing Recovery For Adaptive Decision Feedback Equalization of The Magnetic Stor- 
age Channel". Glot>ecom'90 IEEE Global Telecommunications Conference 1990, San Diego. CA. Nov. 1990, pp. 1794- 
1799; and Abbott et al. "Performance of Digital Magnetic Recording with Equalization and Offtrack Interference", IEEE 
Transactions on Magnetics, Vol. 27, No. 1 , Jan. 1991 ; and Cioffi et al, "Adaptive Equalization in Magnetic-Disk Storage 
Channels", IEEE Communication Magazine, Feb. 1990; and Roger Wood, "Enhanced Decision Feecback Equalize- 
rs tion", Inter mag'90. 

Similar to conventional peak detection systems, sampled amplitude detection requires timing recovery in order to 
correctly extract the digital sequence. Rather than process the continuous signal to align peaks to the center of bit cell 
periods as in peak detection systems, sampled amplitude systems synchronize the pulse samples to the baud rate. In 
conventional sampled amplitude read channels, timing recovery synchronizes a sampling clock by minimizing an error 
30 between the signal sample values and estimated sample values. A pulse detector or slicer determines the estimated 
sample values from the read signal samples. Even in the presence of ISI the sample values can be estimated and, 
together with the signal sample values, used to synchronize the sampling of the analog pulses in a decision-directed 
feedback system. 

A phase-locked-loop (PLL) normally implements the timing recovery decision-directed feedback system. The PLL 
35 comprises a phase detector for generating a phase en-or based on the difference between the estimated samples and 
the read signal samples. A PLL loop filter filters the phase error, and the filtered phase error operates to synchronize 
the channel samples to the baud rate. Conventionally, the phase error adjusts the frequency of a sampling clock which 
is typically the output of a variable frequency oscillator (VFO). The output of the VFO controls a sampling device, such 
as an analog-to-digital (A/D) converter, to synchronize the sampling to the baud rate. 
40 As mentioned above, sampled amplitude read channels also commonly employ a discrete time equalizer filter to 
equalize the sample values into a desired partial response (PR4, EPR4, EEPR4. etc.) before sequence detection. To 
this end. adaptive algorithms have been applied to compensate in real time for parameter variations In the recording 
system and across the disk radius. For example, U.S. Pat. No. 5,381 ,359 entitled "Adaptation and Training of Digital 
Finite Impulse Response Filter Within PRML Sampling Data Detection Channel", discloses an adaptive equalizer filter 
45 that operates according to a well known least mean square (LMS) algorithm. 

Wk+1 = Wk-*i*ek*XK. 

where W{^ represents a vector of filter coefficients; ^ is a programmable gain; ei^ represents a sample error between the 
50 filter's actual output and a desired output; and represents a vector of sample values from the filter input. In other 
words, the LMS adaptive equalizer filter is a closed loop feedback system that attempts to minimize the mean squared 
error between an actual output of the filter and a desired output by continuously adjusting the filter's coefficients to 
achieve an optimum frequency response. 

A problem associated with adaptive equalizer filters in sampled amplitude read channels is that the timing recovery 
55 and gain control loops can Interfere with the adaptive feedback loop, thereby preventing the adaptive equalizer filter 
from converging to an optimal state. This non-convergence is manifested by the filter's phase and gain response drifting 
as it competes with the timing and gain control loops. An article by J.D. Coker et al. entitled "Implementation of PRML 
in a Rigid Disk Drive", published in IEEE Transactions on Magnetics, vol. 27. No. 6, Nov. 1991, suggests a three tap 
transversal filter comprising a fixed center tap and symmetric side taps In order to constrain the phase response of the 
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equalizer filter except in terms of a fixed group delay. Constraining the phase response of the adaptive equalizer in this 
manner, however, is a very sub-optimal method for attenuating interference from the timing recovery and gain control 
loops. Furthermorei it significantly reduces control over the adaptive filter's phase response, thereby placing the burden 
of phase compensation on the analog equalizer. 

5 Yet another problem associated with conventional adaptive equalizer filters is an inherent limitation on its order (i.e.. 
the number of coefficients}: because the adaptive equalizer is inside the timing recovery feedback loop, its order must 
be limited to minimize the associated transport delay. Compensating for the deficiencies of the discrete time equalizer 
requires filtering the analog read signal with a higher order analog equalizer prior to the timing recovery loop, which is 
undesirable. 

10 There is, therefore, a need for an adaptive, discrete time equalizer filter in a sampled amplitude read channel hav- 
ing an improved method for constraining the phase and gain response in order to minimize interference from the timing 
recovery and gain control loops. A further aspect of the present invention is to remove the adaptive equalizer, and its 
associated latency, from the timing recovery loop, tiiereby allowing a higher order discrete time filter and a simplified 
analog filter. 

15 

SUMMARY OF THE INVENTION 

A sampled amplitude read channel for magnetic disk recording which asynchronously samples the analog read sig- 
nal, adaptively equalizes the resulting discrete time sample values according to a target partial response, extracts syn- 

20 chronous sample values through interpolated timing recovery, and detects digital data from the synchronous sample 
values using a Viterbi sequence detector is disclosed. To minimize interference from the timing and gain control loc^^s, 
the phase and magnitude response of the adaptive equalizer filter are constrained at a predetermined frequency using 
an optimal orthogonal projection operation as a modification to a least mean square (LMS) adaptation algorithm. Fur- 
ther, with interpolated timing recovery, the equalizer filter and its associated latency are removed from the timing reoov- 

25 ery loop, thereby allowing a higher order discrete time filter and a lower order analog filter. 

BRIEF DESCRIPTION OF THE DRAWINGS 

The above and other aspects and advantages of the present invention will be better understood by reading the fol- 
30 towing detailed description of the invention in conjunction with the drawings, wherein: 

FIG. 1 is a block diagram of a conventional sampled amplitude recording channel. 

FIG. 2A shows an exemplary data format of a magnetic disk having a plurality of concentric tracks comprised of a 
plurality of user data sectors and embedded servo data sectors. 
35 FIG. 26 shows an exemplary format of a user data sector. 

FIG. 3 is a block diagram of the improved sampled amplitude read channel of the present invention comprising 
interpolated timing recovery for generating interpolated sample values and a synchronous data clock for clocking 
operation of a discrete time sequence detector and peak detector. 

FIG. 4A Is a detailed block diagram of the prior art sampling timing recovery comprising a sampling VFO. 
40 FIG. 4B is a detailed block diagram of the interpolating timing recovery of the present invention comprising an inter- 
polator. 

FIG. 5 illustrates the channels samples in relation to the interpolated baud rate samples for the acquisition pream- 
ble. 

FIG. 6 shows an FIR filter implementation for the timing recovery interpolator. 
45 FIG. 7 depicts a cost reduced tnplementation for the timing recovery interpolator. 

FIG. 8A is a block diagram of a conventional adaptive, discrete time equalizer filter in a sampled amplitude read 
channel. 

FIG. 88 shows the adaptive, discrete time equalizer of the present invention. 

FIG. 8C shows an alternative embodiment for the adaptive, discrete time equalizer of tiie present Invention. 
50 FIG. 9A illustrates the present invention adaptive filter's gain response constrained at a normalized frequency of 
1/4T 

FIG. 9B shows the present invention adaptive filter's phase response constrained at a normalized frequency of 
1/4T 

FIG. 1 0 illustrates operation of an orthogonal projection operation of the present invention for constraining the gain 
55 and phase response of the adaptive filter. 

FIG. 1 1 A shows an implementation for a reduced cost orthogonal projection operation. 
FIG. 1 1 B shows an alternative embodiment for the reduced cost orthogonal projection operation. 
FIG. lie illustrates an implementation for a gradient averaging circuit used in the reduced cost orthogonal projec- 
tion operation of FIG. 1 1 B. 
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FIG. 1 1 D shows yet another alternative embodiment for the reduced cost orthogonal projection operation of the 
present invention. 

DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENT 

5 

Conventional Sampled Amplitude Read Channel 

Referring now to FIG. 1 , shown is a detailed block diagram of a conventional sampled amplitude read channel. Dur- 
ing a write operation, either user data 2 or prean^le data from a data generator 4 (for example 2T preamble data) is 

10 written onto the media. An RLL encoder 6 encodes the user data 2 into a binary sequence b(n) 8 according to an RLL 
constraint. A precoder 10 precedes the binary sequence b(n) 8 in order to compensate for the transfer function of the 
recording channel 18 and equalizer filters to form a preceded sequence ~b(n) 12. The preceded sequence -b(n) 12 is 
converted into symbols a(n) 16 by translating 14 '-b(N) = 0 into a(N) = -1. and ~b(N) = 1 into a{N) = +1 . Write circuitry 
9, responsive to the symbols a(n) 16, modulates the current in the recording head coil at the baud rate 1/T to record the 

15 binary sequence onto the media. A frequency synthesizer 52 provides a baud rate write clock 54 to the write circuitry 9 
and is adjusted by a channel data rate signal (CDR) 30 according to the zone the recording head is over. 

When reading the recorded binary sequence from the media, timing recovery 28 first lod^ to the write frequency 
by selecting, as the input to the read channel, the write clock 54 through a multiplexor 60. Once locked to the write fre- 
quency, the multiplexor 60 selects the signal 19 from the read head as the input to the read channel in order to acquire 

20 an acquisition preamble recorded on the disk prior to the recorded user data. A variable gain amplifier 22 adjusts the 
amplitude of the analog read signal 58, and an analog filter 20 provides initial equalization toward the desired response 
as well as attenuating aliasing noise. A sampling device 24 samples the analog read signal 62 from the analog filter 20, 
and a discrete time equalizer filter 26 provides further equalization of the sample values 25 toward the desired 
response. In partial response recording, for example, the desired response is often selected from Table 1 . 

25 The discrete equalizer filter 26 may be implemented as a real-time adaptive filter which compensates for parameter 
variations over the disk radius (i.e., zones), disk angle, and environmental conditions such as temperature drift. To this 
end, the filter 26 receives estimated sample values B143 generated by the timing recovery circuit 28; the estimated 
samples being input into an adaptive feedback loop and used to generate sample errors. The adaptive feedback loop 
conventionally employs a least mean square (LMS) algorithm to adapt the filter coefficients (i.e., it adapts the frequency 

30 and phase response of the filter) until a minimum sample error is achieved. Operation of a conventional adaptive equal- 
izer filter is discussed In greater detail below. 

After equalization, the equalized sample values 32 are applied to a decision directed gain control 50 and timing 
recovery 28 circuit for adjusting the amplitude of the read signal 58 and the frequency and phase of the sampling device 
24, respectively. Timing recovery adjusts the frequency of sampling device 24 over line 23 in order to synchronize the 

35 equalized samples 32 to the baud rate. Frequency synthesizer 52 provides a course center frequency setting to the tim- 
ing recovery circuit 28 over line 64 in order to center the timing recovery frequency over temperature, voltage, and proc- 
ess variations. The channel data rate (CDR) 30 signal adjusts a frequency range of the synthesizer 52 according to the 
data rate for the current zone. Gain control 50 adjusts the gain of variable gain amplifier 22 over line 21 in order to match 
the magnitude of the channel's frequency response to the desired partial response. 

40 The equalized samples Y(n) 32 are also sent to a discrete time sequence detector 34. such as a maximum likeli- 
hood (ML) Viterbi sequence detector, which detects an estimated binary sequence ^(n) 33 from the sanple values. An 
RLL decoder 36 decodes the estimated binary sequence ^b(n) 33 from the sequence detector 34 into estimated user 
data 37. A data sync detector 66 detects the sync mark 70 (shown in FIG. 2B) in the data sector 15 in order to frame 
operation of the RLL decoder 36. In the absence of errors, the estimated binary sequence ^b(n) 33 matches the 

45 recorded binary sequence b(n) 8, and the decoded user data 37 matches the recorded user data 2. 

Data Format 

FIG. 2A shows an ^cemplary data format of a magnetic media comprising a series of concentric data tracks 13 
50 wherein each data track 13 comprises a plurality of sectors 15 with embedded servo wedges 17. A servo controller (not 
shown) processes the servo data in the servo wedges 17 and, in response thereto, positions the read/write head over 
a desired track. Additionally, the servo controller processes servo bursts within the servo wedges 17 to keep the head 
aligned over a centerline of the desired track while writing and reading data. The servo wedges 17 may be detected by 
a simple discrete time pulse detector or by the discrete time sequence detector 34. If the sequence detector 34 detects 
55 the servo data, then the format of the servo wedges 1 7 includes a preamble and a sync mark, similar to the user data 
sectors 15. 

FIG. 2B shows the format of a user data sector 15 comprising an acquisition preamble 68, a sync mark 70, and 
user data 72. Timing recovery uses the acquisition preamble 68 to acquire the correct sampling frequency and phase 
before reading the user data 72, and the sync mark 70 demarks the beginning of the user data 72 (see co-pending U.S. 
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patent application 08/313.491 entitled "Improved Timing Recovery For Synchronous Partial Response Recording"). 

To increase the overall storage density, the disk is partitioned into an outer zone 1 1 comprising fourteen data sec- 
tors per track, and an inner zone 27 comprising seven data sectors per track. In practice, the disk Is actually partitioned 
into several zones with a different number of sectors in each zone, and the data recorded and detected at a different 
5 data rate in each zone. 

Improved Sampled Amplitude Read Channel 

FIG. 3 shows the improved sampled amplitude read channel of the present invention wherein the conventional 
10 sanpled timing recovery 28 of FIG. 1 has been replaced by interpolated timing recovery B100. In addition, the write 
frequency synthesizer 52 generates a baud rate write clock 54 applied to the write circuitry 9, and an asynchronous 
read clock 54 for clocking the sampling device 24, the discrete time equalizer filter 26, and the interpolated timing recov- 
ery B100 at a frequency relative to the current zone (CDR 30). In an alternative embodiment, a first frequency synthe- 
sizer generates the write clock, and a second frequency synthesizer generates the read clock. 
15 The discrete equalizer filter B103 is real-time adaptive, receiving Interpolated sample values 81 02 and estimated 
sample values B143 from the interpolated timing recovery circuit B100 for use in a modified least mean square (LMS) 
algorithm which constrains the filter's gain and phase response according to the present invention, the details of which 
are set forth below. 

The Interpolated timing recovery B100 interpolates the equalized sample values 32 to generate interpolated sam- 
20 pie values B102 substantially synchronized to the data rate of the current zone. A discrete time sequence detector 34 
detects an estimated binary sequence 33 representing the user data from the interpolated sample values B102. The 
interpolated timing recovery B100 circuit generates a synchronous data clock B104 for clocking operation of the gain 
control 50, discrete time sequence detector 34, sync mark detector 66 and RLL decoder 36. 

25 Conventional Timing Recovery 

An overview of the conventional sampling timing recovery 28 of FIG. 1 is shown in FIG. 4A. The output 23 of a var- 
iable frequency oscillator (VFO) B1 64 controls the sampling clock of a sampling device 24 which is typically an analog- 
to-digital converter (A/D) In digital read channels. A multiplexor B159 selects the unequallzed sample values 25 during 

30 acquisition and the equalized sample values 32 during tracking, thereby removing the discrete equalizer filter 26 from 
the timing loop during acquisition In order to avoid its associated latency. A phase error detector B155 generates a 
phase error in response to the sample values received over line B149 and estimated sample values -Y^ from a sample 
value estimator B141 , such as a slicer in a d=0 PR4 read channel, over line B143. A loop filter B160 filters the phase 
error to generate a frequency offset A/ B167 that settles to a value proportional to a frequency difference between ttie 

35 sampling clock 23 and the baud rate. The frequency offset Af B1 67, together with the center frequency control signal 
64 from the frequency synthesizer 52. adjust tiie sampling clock 23 at the output of the VFO B164 in order to synchro- 
nize the sampling to the baud rate. 

A zero phase start B1 62 circuit suspends operation of tiie VFO B1 64 at the beginning of aoquisition in order to min- 
imize the initial phase error between the sampling clock 23 and the read signal 62. This is achieved by disabling the 

40 VFO B164, detecting a zero crossing in the analog read signal 62, and re-enabling the VF0 164 after a predetermined 
delay between the detected zero crossing and tiie first baud rate sample. 

The estimated sample values B143 at The output of the slicer B141 are also input into the discrete time equalizer 
fitter 26 of FIG. 1 for use In a conventional least mean square (LMS) adaptation algorithm as is described in more detail 
below. 

45 

Interpolated Timing Recovery 

The interpolated timing recovery B100 of the present Invention is shown in FIG. 4B. The VFO B164 in the conven- 
tional timing recovery of FIG. 4A is replaced with a modulo-Ts accumulator B120 and an interpolator B122. In addition, 

50 an expected sample value generator B151 . responsive to Interpolated sample values B102, generates expected sam- 
ples used by the phase error detector B155 to compute the phase error during acquisition. A multiplexor B153 
selects the estimated sample values -Yi^^^ from the slicer B141 for use by the phase error detector B155 during track- 
ing. The data clock B104 Is generated at the output of an AND gate B126 In response to the sampling clock 54 and a 
mask signal B1 24 from the modu!o-Ts accumulator B1 20 as discussed in further detail below. The phase error detector 

55 B155 and the slicer B141 process interpolated sample values B102 at the output of the interpolator B122 rather than 
the channel sample values 32 at the output of the discrete equalizer filter 26 as in FIG. 4A. A PID loop filter B161 con- 
trols the closed loop frequency response similar to the loop filter B1 60 of FIG. 4A. 

The interpolated sample values Y^^^ B102 and the estimated sample values -V^^^ from the slicer B141 are Input 
into the adaptive, discrete equalizer filter B103 of FIG. 3 for use by a modified least mean square (LMS) algorithm, tiie 
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details of which are set forth below. 

In the interpolated timing recovery of the present invention, locking a VFO to a reference frequency before acquir- 
ing the preamble is no longer necessary; multiplexing 60 the write clock 54 into the analog receive filter 20 (as in FIG. 
1) is not necessary. Further, the sampling device 24 and the discrete equalizer filter 26, together with their associated 

5 delays, have been removed from the timing recovery loop; it is not necessary to multiplex B159 around the equalizer 
filter 26 between acquisition and tracking. However, it is still necessary to acquire a preamble 68 before tracking the 
user data 72. To this end, a zero phase start circuit B163 minimizes the Initial phase error between the interpolated 
sample values and the baud rate at the beginning of acquisition similar to the zero phase start circuit B162 of FIG. 4A. 
However, rather than suspend operation of a sampling VFO B1 64, the zero phase start circuit B1 63 for interpolated tim- 

w ing recovery computes an initial phase error t from the A/D 24 sample values 25 and loads this initial phase error into 
the modulo-Ts accumulator B1 20. 

For more details concerning the PID loop filter B1 61 , phase error detector B1 55. expected sample generator B1 51 . 
and slicer B141 , refer to the above referenced co-pending U.S. patent applications "Sampled Amplitude Read Channel 
Comprising Sample Estimation Equalization, Defect Scanning, Channel Quality, Digital Servo Demodulation, PID Filter 

IS for Timing Recovery, and DC Offset Control" and "Improved Timing Recovery For Synchronous Partial Response 
Recording." A detailed description of the modulo-Ts accumulator B120, data clock B104, and interpolator B122 is pro- 
vided in the following discussion. 

Interpolator 

20 

The interpolator B122 of FIG. 48 is understood with reference to FIG. 5 which shows a sampled 2T acquisition pre- 
amble signal B200. The target synchronous sample values B102 are shown as black circles and the asynchronous 
channel sample values 32 as vertical arrows. Beneath the sampled preamble signal is a timing diagram depicting the 
corresponding timing signals for the sanpling clock 54, the data clock B1 04 and the mask signal B1 24. As can be seen 
25 in FIG. 5, the preamble signal B200 is sampled slightly faster than the baud rate (the rate of the target values). 

The function of the interpolator is to estimate the target sample value by interpolating the channel sample values. 
For illustrative purposes, consider a simple estimation algorithm, linear interpolation: 

Y(N-1) = x(N-1) + T • (x(N) - x(N-1)); where: (1) 

30 

x(N-1) and x(N) are the channel samples surrounding the target sample; and x is an interpolation interval proportional 
to a time difference between the channel sample value x(N-1) and the target sample value. The interpolation interval t 
is generated at the output of modulo-Ts accumulator B120 which accumulates the frequency offset signal B167 at 
the output of the PID loop filter 81 61 : 

35 

X = (ZAO MOD TS; where: (2) 

Ts is the sampling period of the sanpling clock 54. Since the sampling clock 54 samples the analog read signal 62 
slightly faster than the baud rate, it is necessary to mask the data clock every time the accumulated frequency offset Af, 

40 integer divided by Ts, increments by 1 . Operation of the data clock B104 and the mask signal B124 generated by the 
modulo-Ts accumulator B120 is understood with reference to the timing diagram of FIG. 5. 

Assuming the interpolator implements the simple linear equation (1) above, then channel sample values B202 and 
B204 are used to generate the interpolated sample value corresponding to target sample value B206. The interpolation 
interval x B208 is generated according to equation (2) above. The next interpolated sample value corresponding to the 

45 next target value B21 0 is computed from channel sanrple values B204 and B21 2. This process continues until the inter- 
polation interval x B214 would be greater than Ts except that it "wraps" around and is actually x B216 (I.e., the accumu- 
lated frequency offset Af, integer divided by Ts, increments by 1 causing the mask signal B124 to activate). At this point, 
the data clock B104 is masked by mask signal B124 so that the interpolated sample value corresponding to the target 
sample value B220 is computed from channel sample values B222 and B224 rather than channel sample values B218 

50 and B222. 

The simple linear interpolation of equation (1) will only work if the analog read signal is sampled at a much higher 
frequency than the baud rate. This is not desirable since operating the channel at higher frequencies increases its com- 
plexity and cost. Therefore, In the preferred embodiment the interpolator B1 22 is implemented as a filter responsive to 
more than two channel samples to compute the interpolated sample value. 
55 The ideal discrete time phase interpolation filter has a flat magnitude response and a constant group delay of x: 

C,(eh = e^"^ (3) 

which has an ideal impulse response: 
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sinc(7c • (n-x/Tg)). (4) 

Unfortunately; the above non-causal Infinite impulse response (4) cannot be realized. Therefore, the impulse response 
of the Interpolation filter is designed to be a best fit approximation of the ideal Impulse response (4). This can be accom- 
5 plished by minimizing a mean squared error between the frequency response of the actual interpolation filter and the 
frequency response of the ideal interpolation filter (3). This approximation can be improved by taking into account the 
spectrum of the input signal, that is. by minimizing the mean squared error between the input spectrum multiplied by 
the actual interpolation spectrum and the input spectrum multiplied by the ideal interpolation spectrum: 

10 C,(e^^)X(eh ■ C^(ehx(e^^); where: (5) 

C^(e^^) is the spectrum of the actual interpolation filter; and X(e^^) is the spectrum of the input signal. From equation 
(5), the mean squared error is represented by: 

15 iz 

=2^J|C,(eh-en'|X(ehl'cffi);where: (6) 



20 X(e^) is the spectrum of the read channel (e.g., PR4» EPR4, EEPR4 of Table 1 or some other partial response spec- 
trum). 

In practice, the above mean squared error equation (6) Is modified by specifying that the spectrum of the input sig- 
nal is bandlimitted to some predetermined constant 0 ^ co ^ an where 0 < a < 1 ; that is: 

25 \X(e^'^)\ = 0. for ^ an. 

Then equation (6) can be expressed as: 



30 



50 



an 



E?.„ = £ J |C,(eh ■ en'm^"')\'d^. (7) 

-an 



The solution to the minimization problem of equation (7) involves expressing the actual interpolation filter in terms of its 
35 coefficients and then solving for the coefficients that minimize the error in a classical mean-square sense. 
The actual interpolation filter can be expressed as the FIR polynomial: 

n=R-1 



2R is the number of taps in each Interpolation filter and the sample period Ts has been normalized to 1 . A mathematical 
derivation for an interpolation filter having an even number of coefficients is provided below. It is within the ability of 
45 those skilled in the art to modify the mathematics to derive an interpolation filter having an odd number of coefficients. 
Substituting equation (8) into equation (7) leads to the desired expression in terms of the coefficients C^(n): 

E'a = i J I Z C,(n)e-""" - e""^|^|X(e"")|^d«>. (9) 



2n 

-aji n=-R 



The next step Is to take the derivatives of equation (9) with respect to the coefficients C^(n) and set them to zero: 

ac'(n ) =Oforno = -R 0.1 R-1. (10) 

After careful fnanv}ulation, equation (10) leads to: 
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an 



J 



•an 



-/-nsR-l \ 

y C^(n)cos(ffl(no-n)) 



• C08( o)(n<, + x))l|X(e''')|^cyo» = 



(11) 



10 



for Hn s -R 0, 1 R-1 . 



Defining ^(r) as: 



an 



Hf) = J IXCe'"")! ^ cos((or) da 



(12) 



•an 



IS 



and substituting equation (12) into equation (1 1) gives: 



n=R-1 

X C,(n)(|>{n-n^) = (t>{n^ + x) 

n=-R 



(13) 



20 



25 



30 



for n^ = -R, .... 0, 1, R-1 . 

Equation (13) defines a set of 2R linear equations in terms of the coefficients C^(n), Equation (13) can be expressed 
more compactly In matrix form: 



is a column vector of the form: 



Oj is a Toeplitz matrix of the form: 



<D-pC<5 = 0^: where: 



C, = [c,(-R) c,{0) c,(R-1)]' 



35 



40 



<|)(0) <|)(1) 



<{)(2R-1 



<|){2R-1) 



45 



50 



55 



and is a column vector of the form: 



= [(K-R + x) <|>(x). +(1 + x), .... cMR-1 + x)]\ 



The solution to equation (14) Is: 



= <Dj "^o^; where: 



(14) 



(15) 



Of^ is an inverse matrix that can be solved using well known methods. 

Table B2 shows example coefficients C^(n) calculated from equation (15) with 2R=6. a=0.8 and X(e^'^)=PR4 . The 
implementation of the six tap FIR filter is shown in FIG. 6. A shift register B250 receives the channel samples 32 at the 
sampling clock rate 54. The filter coefficients C^(n) are stored in a coefficient register file B252 and applied to corre- 
sponding multipliers according to the current value of x 8128. The coefficients are multiplied by the channel samples 
32 stored in the shift register 8250. The resulting products are summed 8254 and the sum stored in a delay register 
B256. The coefficient register file B252 and the delay register 8256 are clocked by the data clock B104 to implement 
the masking function described above. 

In an alternative embodiment not shown, a plurality of static FIR filters, having coefficients that conrespond to the 
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different values of t. filter tlie sample values In the shift register B250. Each filter outputs an interpolation value, and the 
current value of the interpolation interval t B128 selects the output of the corresponding filter as the output B102 of the 
interpolator 81 22. Since the coefficients of one filter are not constantly updated as in FIG. 6. this multiple filter embod- 
iment increases the speed of the interpolator B122 and the overall throughput of the read channel. 

Cost Reduced Interpolator 

Rather than store all of the coefficients of the interpolation filters in memory, in a more efficient, cost reduced imple- 
mentation the coefficient register file B252 of FIG. 6 computes the filter coefficients C^(n) in real time as a function of x. 
For example, the filter coefficients C^(n) can be computed in real time according to a predetermined polynomial in x 
(see. for example. U.S. Pat. No. 4,866,647 issued to Farrow entitled, "A Continuously Variable Digital Delay Circuit," the 
disclosure of which is hereby incorporated by reference). An alternative, preferred embodiment for computing the filter 
coefficients in real time estimates the filter coefficients according to a reduced rank matrix r^resentation of the coeffi- 
cients. 

The bank of filter coefficients stored in the coefficient register file B252 can be represented as an MxN matrix Pi^xN* 
where N is the depth of the interpolation fitter (i.e., the number of coefficients C^(n} in the inpulse response computed 
according to equation (15)) and M is the number of interpolation intervals (i.e., the number of x intervals). Rather than 
store the entire A^xn matrix in memory a more efficient, cost reduced implementation is attained through factorization 
and singular value decomposition (SVD) of the A^xn nnatrix. 

Consider that the A^^n matrix can be factored into an Fj^^n and GnxN matrix, 

Then a reduced rank approximation of the Amxn matrix can be formed by reducing the size of the F^xn and GnxN matri- 
ces by replacing N with L where L<N and, preferably, L«N. Stated differently, find the F^xl and Glxm matrices whose 
product best approximates the A|^xn matrix, 

^MxN * ^MxL * Qlxn ' 

The convolution process of the interpolation filter can then be carried out, as shown in FIG. 7, by implementing the Glxn 
matrix as a bank of FIR filters B260 connected to receive the channel sample values 32, and the Fiv^xl matrix imple- 
mented as a lookup table B262 indexed by x B128 (as will become more apparent in the following discussion). Those 
skilled in the art will recognize that, in an alternative embodiment, the A^xn matrix can be factored into more than two 
matrices (i.e., A « FGH...). 

The preferred method for finding the F^xl and G|_xn matrices is to minimize the following sum of squared en^ors: 

M N 

S I(Ajn-(NxL • Glxn)>i)' (16) 

j=1n*1 



The solution to equation (16) can be derived through a singular value decomposition of the A^xn matrix, comprising the 
steps of: 

1 . performing an SVD on the A|^xN matrix which gives the following unique factorization (assuming M^N): 

Amxn = U MXN • D NXN • Vnxn where: 

Umxn a MxN unitary matrix; 

Dnxn is a NxN diagonal matrix {ai, a^} where aj are the singular values of A^^xn» and <ti ^ 02... ^ afg ^ 
0; and 

Vnxn is a NxN unitary matrix; 

2. selecting a predetermined L number of the largest singular values a to generate a reduced size diagonal matrix 
Dlxl: 



10 



EP0805447 A2 



5 



□lxl =Diag {aj.Cj, —tC^} = 



10 



Oj 0 ... 0 

0 02 0... : 

: ... . 0 

0 ... 0 Cj 



IS 



20 



2S 



3. extracting the first L columns from the U^xN i^atrix to form a reduced U^xL matrix: 



U 



MxL - 



^1,1 ••• ^1,L 



• • • • « 



• • • • • • 



U 



1,N 



u 



M,N 



30 



35 



40 



4. extracting the first L rows from the V(sjxim matrix to form a reduced V|_xn matrix: 



V 



LxN 




45 

5. defining the F^xL ^"d Glxn matrices such that: 

50 (for example, let F^^^l = ^mxl * ^lxl ^^<^ Glxn = Vlxn)- 

In the above cost reduced polynomial and reduced rank matrix embodiments, the interpolation filter coefficients 
C^(n) are computed in real time as a function of x\ that is. the filter's impulse response h(n) is approximated according to: 

55 L 

h(n.T) = c,(n) = X ' K'.^): where: (17) 

i« 1 
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f(i.t) is a predetermined function in x (e.g., polynomial in x or x indexes the above Fi^^l Tiatrix); L is a degree which 
determines the accuracy of the approximation (e.g., the order of the polynomial or the column size of the above F^xl 
matrix); and Qi(n) is a predetermined matrix (e.g., the coefficients of the polynomial or the above Gun nfiatrix). As L 
increases, the approximated filter coefficients C^{n) of equation (17) tend toward the ideal coefficients derived from 
5 equation (1 5). It follows from equation (1 7) that the output of the interpolation filter Y(x) can be represented as: 

N L 

Y(x)= U(x-n) 5^G,(n) • f(i.x) (18) 

10 

where U(x) are the channel sample values 32 and N is the number of interpolation filter coefficients C^(n). 

Referring again to FIG. 6, the coefficient register file can compute the interpolation filter coefficients C^(n) according 
to equation (1 7) and then convolve the coefficients C^(n) with the channel samples U(x) 32 to generate the interpolated 
sample values B102 synchronized to the baud rate. However, a more efficient implementation of the interpolation filter 
15 can be achieved by rearranging equation (18): 

L N 

. Y{x)= X^''-^) SGiW • U(x-n) (19) 

1=1 n=1 

20 

FIG. 7 shows the preferred embodiment of the interpolation filter according to equation (19). In the polynomial 
embodiment, the function of x is a polynomial in x, and the matrix Gj(n) are the coefficients of the polynomial. And in the 
reduced rank matrix embodiment, the function of x is to index the above F|vixl matrix B262. and the second summation 
25 In equation (19), 

N 

XQ|(n) • U(x-n) 

n= 1 

30 

is implemented as a bank of FIR filters B260 as shown in FIG. 7. Again, in equation (19) L is the depth of the approxi- 
mation function f(i,x) (e.g., the order of the polynomial or the column size of the above F^^l matrix) and N is the depth 
of the interpolation filter's impulse response (i.e., the number of coefficients in the impulse response). It has been deter- 
35 mined that N=8 and Ls3 provides the best performance/cost balance; however, these values may increase as IC tech- 
nology progresses and the cost per gate decreases. 

Conventional Adaptive Equalizer 

40 FIG. 8A illustrates a prior art adaptive, discrete time equalizer that operates according to the well known least mean 
square (LMS) algorithm, 

45 or alternatively, 

Wk+1 = Wj^-^-X^-eK 

where represents a vector of FIR filter coefficients; is a programmable gain; ei^ represents a sample error (or vec- 
50 tor of sample errors e^) between the FIR filter's actual output and a desired output; and represents a vector of sam- 
ples values (or a scalar from the FIR filter input. To better understand operation of the present invention, the second 
representation of the LMS algorithm is used throughout this disclosure. 

The desired filter output is the estimated sample values -Y^ at the output of slicer 8141 . The estimated sample val- 
ues -\ are subtracted from the FIR filter's output Y^ to generate the sample error ek. The LMS algorithm attempts to 
55 minimize the sample error in a least mean square sense by adapting the FIR filter coefficients; that is, it adjusts the FIR 
filter's gain and phase response so that the overall channel response adequately tracks the desired partial response 
(e.g., PR4, EPR4. EEPR4, etc.). 

As previously mentioned, interference from the timing recovery 28 and gain control 50 loops can prevent the adap- 
tive, discrete time equalizer 26 from converging to an optimal state. For example, a phase adjustment in the adaptive 
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filter 26 can affect the sampling phase error for timing recovery 28. Timing recovery 28 compensates for the filter's 
phase adjustment by adjusting its sampling phase; this adjustment can result in yet another phase adjustment by the 
adaptive equalizer 26. Thus, the phase response of the adaptive equalizer may never converge. Similarly, the gain con- 
trol loop 50 can interfere with the gain response of the adaptive filter 26 and prevent it from converging. 

5 

Constrained Adaptive Equalizer 

FIG. 8B illustrates operation of the adaptive, discrete time equalizer filter B103 of the present invention. Sample val- 
ues from the A/D 24 are input over line 25 to a discrete time FIR filter CI 00 comprising a predetermined number of ooef- 
10 ficients, the values of which determine the filter's gain and phase response. Because the FIR filter CI 00 operates on 
the sample values prior to the interpolated timing recovery loop 81 00, its order can be increased over the prior art with- 
out adversely affecting the latency of timing recovery (i.e., the number of filter coefficients can be increased). 

The output Y|( 32 of the FIR filter C100 is input into the interpolator 8122 for generating the interpolated sample 
values 8102. The interpolated sample values Y^^.^ 8102 are input into a slicer 81 41 (FIG. 4B) which generates 
15 estimated sample values -Y^+^j. The estimated sample values ~\+^ are subtracted from the interpolated sample val- 
ues at adder CI 02 to generate a sample error value e^+x CI 04 that is synchronized to the baud rate rather than 
the sample rate. Because the LMS algorithm operates on sample values X|^ at the sample rate, it is necessary to con- 
vert the error value 6^+^ CI 04 into an error value ei^ C112 synchronous to the sample rate. This is accomplished by an 
interpolation circuit C106 which computes an interpolated error value C112 from the baud rate error values e^+t 
20 CI 04. Preferably, the error value interpolation circuit CI 06 is implemented as a first order linear interpolation, but it may 
be a sinple zero order hold, or a more complex Interpolation filter as describe above. 

FIG. 80 shows an alternative embodiment for generating the error value e^^. As illustrated, the estimated ideal sam- 
ple values -Y}c^ from the slicer 8141 are interpolated by the interpolation circuit CI 06 to generate estimated ideal sam- 
ple values -Y}^ CI 20 which are subtracted CI 02 from the equalized sample values Y^ 32 at the output of the FIR filter 
25 CI 00 to generate the error value e|^ 

In both embodiments, the error value ei^ C1 12 is input into a modified LMS circuit C1 14 which computes updated 
filter coefficients W|^ C1 16 according to, 

WK.i = W,-^-Pv,V2^.(XK-eK) 

30 

where the operation Pv-i vg-*- is an orthogonal projection operation which constrains the gain and phase response of the 
FIR filter C100 in order to attenuate interference from the gain and timing loops. 

Operation of the orthogonal projection operation Pv^vg-^ will now be described in relation to the gain and phase 
response of the FIR filter CI 00. FIG.s 9A and 9B show an example gain and phase response, respectively, for the FIR 

35 filter CI 00 of the present invention. The gain and phase response vary over time as the filter adapts to parameter 
changes in the recording system; that is, the filter continuously adjusts the channel's overall frequency response so that 
it matches the desired partial response (PR4, EPR4, EEPR4, etc.) as best possible. In the present invention. Interfer- 
ence from the timing recovery and gain control loops is attenuated by constraining the gain and phase response of the 
FIR filter CI 00 at a predetermined frequency using an orthogonal projection operation Pv^vg-^. 

40 Referring again to FIG. 9A, the gain (magnitude) response of the FIR filter CI 00 has been constrained to a prede- 
termined value (denoted by g) at the normalized frequency of 0.5 (1/4Ts). Sirnilarly, the phase response of the FIR filter 
CI 00 has been constrained to kn at the normalized frequency of 0.5 as shown in FIG. 9B. In effect, these constraints 
allow the gain and phase to vary (adapt) at all frequencies except at the normalized frequency of 0.5, thereby constrain- 
ing the filter's frequency response in a manner that attenuates interference from the gain and timing loops. The gain 

45 constraint g is relatively arbitrary except that it is selected to optimize the dynamic range of the filter's coefficients. How- 
ever, constraining the filter's response at the normalized frequency of 0.5 and selecting a phase constraint of kn 
reduces the complexity of the orthogonal projection operation PV1V2-'", and simplifies implementation of the zero phase 
start circuit 8163 (FIG. 4B) and sync mark detector 66 (FIG. 1). 

As mentioned above, the zero phase start circuit 8163 (FIG. 4B) minimizes the initial phase error between the inter- 

50 polated sample values and the baud rate at the beginning of acquisition by computing an initial phase error t from the 
A/D 24 sample values 25 and then loading this initial phase error into the modulo-Ts accumulator 8120. To compute 
the Initial phase error x, the zero phase start circuit 8163 must take into account the phase delay of the adaptive equal- 
izer filter 8103 since interpolated timing recovery 81 00 operates on the equalized samples 32. not the A/D 24 samples 
25. With the acquisition preamble 68 (FIG. 2B) having a frequency of 1/4T (i.e.. 0.5 normalized), constraining the phase 

55 response of the adaptive equalizer 81 03 at the preamble frequency (1/4Ts) fixes the phase delay of the equalizer 8103 
during acquisition, thereby allowing the zero phase start circuit 8163 to accurately compute the initial phase error x. 
Furthermore, since the phase constraint at the preamble frequency is fixed at kn, the phase delay of the equalizer 8103 
will either be zero or 180° (i.e., the adjustment to the initial phase error is nothing or a sign change). 

Constraining the phase response of the adaptive equalizer 8103 to kn at the preamble frequency also simplifies 
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10 



implementation of the sync mark detector 66 (FIG. 1) in sampled amplitude read channels that use the conventional 
synchronous sampling timing recovery 28 (FIG. 4A) rather than Interpolated timing recovery B100 (FIG. 4B). Operation 
of the sync mark detector 66 is described in detail in the at>ove referenced co-pending U.S. Patent Application 
08/533,797 entitled "Improved Fault Tolerant Sync Mark Detector For sampled Amplitude Magnetic Recording." As 
described therein, the sync mark detector 66 is enabled coherent with the end of the acquisition preamble 68 and rela- 
tive to the transport delay from the output of the A/D 24 to the sync mark detector 66. With an adaptive equalizer, the 
transport delay will vary unless the filter's phase response is constrained at the acquisition preamble frequency by. for 
example, using an orthogonal projection operation Pyf^v2^ of the present invention. 

Turning now to the implementation details of the orthogonal projection operation Pv^V2-^. the equalizer's frequency 
response is 



15 



C(eJ'"') = SCKe 

k 



-jk27ilT 



where are the coefficients of the equalizer's impulse response. At the preamble frequency (1/4T), the equalizer's fre- 
quency response is 



20 



C(e') = XCKe ' 



25 where the sampling period has been normalized to T=1 . In matrix form, the equalizer's frequency response at the pre- 
amble frequency is, 



30 



35 



40 



f 
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e 



^2 



e 



-I 



e 



(j)-' 




1 

-J 
-1 


• 

.(j) 




• 

J 

■ 



45 



Those skilled in the art will recognize that shifting the time base will lead to four different, but functionally equivalent, 
frequency responses at the preamble frequency (i.e.. [l.-j.-l j.1,...]C, [-1,j,1.-j,...]C and D.1.i-1. -]S)- Con- 

straining the phase response of the equalizer BIOS to an integer multiple of tc at the preamble frequency (1/4T) inrpiies 
that the imaginary component of its frequency response is zero, 



50 



0 
-1 
0 

I 



55 



If the imaginary component of the frequency response is constrained to zero, as described above, then constraining the 
magnitude of the equalizer to g at the preamble frequency (1/4T) implies that the real component of the frequency 
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T 



Therefore, the equalizer's coefficients C|( must be constrained to satisfy the following two conditions: 

c'^-Vi =0andC"^-V2 = g. 

15 

The above constraints are achieved by multiplying the computed gradient X|( • by an orthogonal projection operation 
Py/^V2^ as part of a modified LMS algorithm C1 14. 

To understand the operation of the orthogonal projection operation, consider an equalizer that comprises only two 
coefficients: Cq and as shown in FIG. 10. The phase constraint condition • =0 implies that the filter coeffi- 

20 cient vector must be orthogonal to . When using an unmodified LMS algorithm to update the filter coefficients, the 
orthogonal constraint is not always satisfi^ as shown in FIG. 10. The present invention, however, constrains the fitter 
coefficients to a subspace <£) which is orthogonal to Vi by multiplying the gradient values X^^o^bya projection oper- 
ation Pv/, where the null space of the projection operation Pv/ is orthogonal to iQ). The updated coefficients corre- 
spond to a point on the orthogonal subspace (Q) closest to the coefficients derived from the unmodified LMS algorithm 

25 as shown in FIG. 10. 

Similar to the phase constraint projection operation Pv^-^. a second orthogonal projection operation Pv2-^ constrains 
the filter coefficients such that the coefficient vector Qj satisfies the above gain constraint: C ^ • V 2 = g . The combined 
orthogonal projection operation Pv^V2'^ eliminates two degrees of freedom in an N-dimensional subspace where N is 
the number of filter coefficients (i.e., the orthogonal projection operation PviV2'^ has a rank of N-2). 
30 An orthogonal projection operation for and V2 can be computed according to 

Pv/ = |.Pv, = l-V,(V,X) V/ (20) 

where Pv^Vg "^=Pv., • Pvg ^ since is orthogonal to V2. The orthogonal projection operation PviV2-'- computed 
35 using the above equation for an equalizer comprising ten filter coefficients is a matrix 



I 
0 
-1 
0 




The above matrix PV1V2''' is an orthogonal projection matrix scaled by 5 (multiplied by 5) so that it contains integer 
valued elements which simplifies multiplying by Xi^ • 6)^ in the LMS update equation, 

55 Wk^i = W,-n-PviV2^-(X,.eK). (21) 



The scaling factor Is taken Into account In the selection of the gain value \i. Constraining the gain to g and the phase to 
kn at the normalized frequency of 0.5 simplifies implementing the orthogonal projection matrix PviV2~'~: half of the ele- 
ments are zero and the other half are either +1,-1, or +4. Thus, multiplying the projection matrix Pv^ V2^ by the gradient 
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' values • requires only shift registers and adders. 

The ACQ/TRK signal shown in FIG. 8B disables adaptation of the FIR filter during acquisition, that is. while acquir- 
ing the acquisition preamble 68 shown in FIG. 2B. Thus, the adaptive equalizer B103 adapts only after acquiring the 
acquisition preamble 68. 

5 

Reduced Cost Orthogonal Constraint Matrix 

Even though the above orthogonal projection matrix Pv^ V2-'- has a simple structure wherein half of the elements are 
zero, it may not be cost effective to directly implement it due to the significant number of shift and accumulate operations 
10 necessary to conpute Pv ^ V2 • (X ^ • e . In order to reduce the cost and complexity, an alternative embodiment of 
the present invention decimates the modified LMS adaptation algorithm C1 14 as illustrated in FIG.s 11 A, 11 B. 11 C and 
11D. 

A mathematical basis for the circuit of FIG. 1 1 A will be discussed before describing the details of operation. Refer- 
ring again to the above equation (20). 

15 

Pv/ = l-Pv, = |.V,(V,X)-V/ 
by combining the above and V2 vectors into a Nx2 matrix 

20 



V = 

25 



(those skilled in the art will recognize that shifting the time base provides four alternatives for the V matrix) then the 
30 operation (Vx^Vx)'^ of equation (20) reduces to 



4' 



Thus, equation (20) reduces to 

Pv-^ = l-1/5-w'^. (22) 

40 

Multiplying both sides of equation (22) by 5 provides 

5'Pv-^ = 5n-w'^. (23) 

45 Referring again to equation (21), 

Wk.i = WK-fA-PviV2^-(Xk-eK) 
setting X • e ,^ = g reduces equation (21) to 

50 

Wk+1 =W^-^-Pv^.g, = W,-^-[5-gK-w^-gJ. (24) 

Defining 5 = v^^g^^ 

55 
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+10-1 0 ... 
0-10+1 ... 



>0 



So 

83 



(25) 



IS 



20 



ss 



that is, 5Q=go-g2+— and 6.,=-g,+g3-... , 
then computing • g^^ 



+1 0 
0 -1 
-1 0 

0 +1 



and computing 5 • g - • g,, provides 



+Sfl 

■So 
+5^ 



(26) 
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(27) 



40 



45 



FIG. 1 1 A implements the above equations (25), (26), (27) and ultimately equation (21) in order to update the coef- 
ficients of the equalizer filter CI 00 according to the modified LMS adaptation algorithm of the present invention. To fur- 
ther reduce the implementation cost, the circuit of FIG. 1 1 A decimates the adaptation algorithm by the number of filter 
coefficients (10 in the example shown); that is, the adaptation algorithm operates only on every tenth sample value 
rather than on every santple value and updates only one filter coefficient per clock period. This Is illustrated by expand- 
ing equation (24) 



50 
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^ (k + 0)0 
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58(k+0)0-S(k + o)0 

58(ic + 3)3'S(t + 3)l 



(28) 



where l<+i is the clock period and g (qj = X j.j • e , . Thus, at each clock period, the next gradient value g(i)j can be com- 
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puted by multiplying a sample value latched every tenth clock cycle by the next error value ej. The new gradient g^^j is 
then used to update the corresponding filter coefficient W(Qj of equation (28). 

Referring now to the circuit shown In FIG. 1 1 A. the sample values X|^ are input into the discrete time equalizer filter 
CI 00 and the equalized sample values Input into the interpolator B122 similar to FIG. 8B. A serial-to-parallel circuit 
C140 converts the interpolated sample values into even and odd subsequences Y2k+x and Y(2k-i)+t. where the 
notation 2k Indicates that two new interpolated sample values are output from the serial-to-parallel circuit C140 at every 
other sample period. A sllcer CI 42 generates the corresponding even and odd estimated subsequences -Y2K+t and 
~Y(2k-i)+T which are subtracted from the Interpolated sample values at respective adders (C144^,C1442) to form an 
even and odd sample error sequences e2kM and e(2k.i)+t. An error value interpolation circuit C146, similar to that of 
FIG. SB, generates the even and odd sample error sequences e2k and e(2k-i} which are synchronized to the A/D sample 
rate. 

As mentioned above, the circuit of FIG. 11 A decimates the adaptation algorithm of the present Invention by the 
number of coefficients In the equalizer filter CI 00. For the example shown, the equalizer filter CI 00 has 1 0 filter coeffi- 
cients; accordingly, a decimate by 10 circuit CI 48 loads a sample value Into a delay register CI 50 every tenth sample 
period. Thereafter, the output of delay register CI 50 is represented by X2k.2].i where: j=1 ->5 incremented by 2 at every 
other sample period. The output of the delay register CI 50 X2k-2j.i is multiplied by the sample errors e2k and e(2k-i) at 
respective multipliers (CI 52^. CI 522) to form the gradient values g2j and g2j.i used in equation (24). 

The gradient values g2j and g2j.i are then shifted into respective shift registers (C154i.C1542). To implement equa- 
tion (27), the gradient values g2j and g2j,i are multiplied by respective alternating ±1 {C156i.C1562) and accumulated 
in respective accumulators (C158^,C1582). ^^ar accumulating 5 gradient values in each accumulator (C158^.C1582). 
the outputs of tiie accumulators (which represent §o and 6^ of equation (25)) are latched by respective decimate by 5 
circuits (C160i,C1602) and the accumulators (C1581.CI582) are reset. The values 5o and 61 are then multiplied by 
respective alternating ±1 (CI 62^, CI 622) to implement equation (26). The gradient values g2| and g2j.i at the outputs of 
the shift register (Cl54i,C1542) are multiplied by 5 (C164^.C1642) and the values 5o and di are subtracted therefrom 
at adders (Cl66i,C1662) in order to Implement equation (27). 

Tofinish the adaptive update algorithm (i.e., to implement equation (28)), tiie output of adders (C166^,C1662) are 
scaled by a gain factor ^ (C16a|,Cl682) which Is reduced by a factor of 5 to account for the scaled up projection oper- 
ator. The output of the gain factor \i (CI 681, CI 682) is subtracted at adders (C170i,C1702) from the corresponding filter 
coefficient {W2|.i,W20 selected by a multiplexor (C172^,C1722) from a bank of registers (Cl74i,C1742). The ADAPTgj 
signal selects the 2^ coefficient from the bank of registers (C174i.C1742) for updating. After subtracting the update 
value, the updated filter coefficient (Cl76i,C1762) is restored to the bank of registers (C174i,Cl742) and used by tiie 
equalizer filter CI 00 during the next clock period to equalize the sample values according to its updated spectrum. 

Decimating the update algorithm by the number of filter coefficients, as in FIG. 1 1 A, decreases the implementation 
complexity but at the cost of slowing convergence of tiie equalizer filter CI 00 toward an optimum setting. This is 
because the decimated update algoritiim does not use all of the sample values to compute the gradients gj. FIGs. 1 1 B 
and 1 10 show an alternative embodiment of the decimated update algorithm of the present invention wherein more of 
the sample values are used to compute the gradients gj, which adds conrplexity but improves performance because tiie 
. equalizer filter CI 00 converges faster. 

The circuit in FIG. 11 B operates similar to the circuit in FIG. 1 1 A described above except for the addition of respec- 
tive gradient averaging circuits (Cl76i,C1762) which compute an averaged gradient g2j.i and g2j over several sample 
values, 

^21^ * N 2 ®2k-n • ^2k-2j-n (29) 
n«0 



where N is a predetermined number of sample values out of tiie number of equalizer filter taps. In one embodiment, all 
of the sample values could be used (I.e., a number of sample values equal to the number filter taps), or in an alternative 
embodiment, a decimated number of sample values could be used in order to reduce the complexity and cost. 

FIG. 1 1C shows an embodiment of the gradient averaging circuits (CI 76^,01 762) wherein the number of sample 
values Is decimated by 2; that is. 5 out of the 10 sample values corresponding to the 10 equalizer filter taps are used to 
compute the averaged gradients g2j.i and g2j. 

In operation, a decimate by 2 circuit CI 78 of FIG. 1 1 B stores every other sample value in the delay register C150. 
The output of the delay register CI 50 is multiplied (CI8O1-I8O5) by the error value e2k and delayed versions of the error 
value e^.^' and tiie results are accumulated in accumulators (C182.f-1825). After accumulating five gradients, the con- 
tents of the accumulators (CI 82^-1825) are transferred to registers (C184^-1845), and the accumulators (C182^-1825) 
are cleared. Then, at every other sample period, the contents of registers (CI 84^-1845) are shifted from left to right (i.e., 
C1845 s CI844: CI 844 s CI843: etc.) and tiie output of register CI845 is the averaged gradient g2j output by the gra- 
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10 



dient averaging circuit (C1 76^, C1 762). The averaged gradients g2j and g2j.i (C186^,C1862) are then used to update the 
coefficients of the equalizer filter C100 using circuitry C188 in the same manner as described with reference to FIG. 
11 A. 

FIG. 1 1 D Illustrates yet another embodiment of the present invention which further reduces the cost and complex- 
ity, as compared to FIG. 1 1 A. by updating the even and odd filter coefficients sequentially. That is, during the first N peri- 
ods (where N is the number of filter coefficients) the circuit of FIG. 11D computes the coefficlerrt update values of 
equation (28) for the even filter coefficients (Wq, W2, W4,..), and then during the next N sample periods it computes the 
update values for the odd filter coefficients (W^. W3. W5,...). Adecimation circuit C190 decimates the error value ej^ by 
two, and approximately half the circuitry as that of FIG. 1 1 A is used to conpute the update values. The decimate by two 
circuit CI 90 Is actually syncopated; that Is. it outputs the error values for ksO, ks2, ks4, k=6. ks8, and then outputs the 
error values for k=9, k=1 1 , k=13, k=15, k=1 7 (assuming the adaptive filter CI 00 comprises ten filter taps). 

Mathematically, operation of the update circuit in FIG. 1 1 D can be described by. 
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35 assuming the adaptive filter CI 00 comprises 1 0 filter taps. 

Obviously, the embodiment of FIG. 11 D will decrease the performance of the adaptive filter ClOO due to the 
decrease in convergence speed. However, the gradient averaging circuit of FIG. 1 1C can be used to improve the per- 
formance of the circuit in FIG. 1 1 D, similar to the embodiment of FIG. 1 1 B. Thus in the enf^xxiiments of FIG. 1 1 A-1 1 D, 
the performance versus cost and complexity varles-the preferred configuration is selected according to the require- 

40 ments of the user. 

Although the interpolated timing recovery and adaptive equalizer of the present invention have been disclosed in 
relation to a d^o PR4 read channel, the principles disclosed herein are equally applicable to other types of sanfipled 
amplitude read channels including d=1 EPR4 or EEPR4 read channels. In a d=1 read channel, for example, the slicer 
B1 41 of FIG. B4A is replaced by a pulse detector as described in the above reference U.S. Patent No. 5.359.63 1 . 
45 Furthermore, the particular constraint frequency of 1/4Ts used in the disclosed embodiment Is not Intended to be 
limiting. Other constraint frequencies could be used without departing from the scope of the present invention. For 
example, a 3T preamble could be used in which case the constraint frequency would be 1/6Ts (if constraining to the 
preamble frequency). 

Still further, those skilled in the art will appreciate the many obvious modifications that are possible to the adaptive 
50 equations disclosed herein. For example, a shift in the time base would^change the FIR filter's magnitude and phase 
response at the constraint frequency which would result in different and V2 matrices. Also, the modified IMS algo- 
rithm, 

Wk*i = W,.^i-PviV2^-(X,.eK) 

56 

could be implemented after rearranging terms, 

WK*i=WK-t'-eK-(PViV2^.X,) 
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or 

WK.i = W,-n.X,.{Pv,v/.eK). 

These, and like modifcations. are within the scope of the present invention. 

The objects of the invention have been fully realized through the embodiments disclosed herein. Those skilled in 
the art will appreciate that the various aspects of the Invention can be achieved through different embodiments without 
departing from the essential function. The particular embodiments disclosed are illustrative and not meant to limit the 
scope of the invention as appropriately construed by the following claims. 



Table 1 



Channel 


Transfer Function 


Dipulse Response 


PR4 


(1-D)(1+D) 


0,1,0. -1,0. 0, 0. ... 


EPR4 


(1-D)(1+D)2 


0. 1. 1, -1, -1.0. 0. ... 


EEPR4 


(1-D) (1+D)3 


0,1.2. 0.-2. -1,0. ... 
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Table B2 



10 



15 



20 



25 



30 



35 



40 



45 



50 



T*32/Ts 


c(-2) 


c(-2) 


c(0) 


C{1) 


c(2) 


c{3) 


0 


0.0000 


-0.0000 


1 .0000 


0.0000 


-0.0000 


0.0000 


1 


0.0090 


-0.0231 


0.9965 


0.0337 


-0.0120 


0.0068 


2 


0.0176 


-0.0445 


0.9901 


0.0690 


-0.0241 


0.0135 


3 


0.0258 


-0.0641 


0.9808 


0.1058 


-0.0364 


0.0202 


4 


0.0335 


-0.0819 


0.9686 


0.1438 


-0.0487 


0.0268 


5 


0.0407 


■0.0979 


0.9536 


0.1829 


-0.0608 


0.0331 


6 


0.0473 


-0.1120 


0.9359 


0.2230 


•0.0728 


0.0393 


7 


0.0533 


-0.1243 


0.9155 


0.2638 


■0.0844 


0.0451 


8 


0.0587 


-0.1348 


0.8926 


0.3052 


-0.0957 


0.0506 


9 


0.0634 


-0.1434 


0.8674 


0.3471 


-0.1063 


0.0556 


10 


0.0674 


-0.1503 


0.8398 


0.3891 


-0.1164 


0.0603 


11 


0.0707 


-0.1555 


0.8101 


0.431 1 


-0.1257 


0.0644 


12 


0.0732 


-0.1589 


0.7784 


0.4730 


-0.1341 


0.0680 


13 


0.0751 


•0.1608 


0.7448 


0.5145 


-0.1415 


0.0710 


14 


0.0761 


-0.1611 


0.7096 


0.5554 


-0.1480 


0.0734 


15 


0.0765 


-0.1598 


0.6728 


0.5956 


-0.1532 


0.0751 


16 


0.0761 


-0.1572 


0.6348 


0.6348 


-0.1572 


0.0761 


17 


0.0751 


-0.1532 


0.5956 


0.6728 


-0.1598 


0.0765 


18 


0.0734 


■0.1480 


0.5554 


0.7096 


-0.1611 


0.0761 


19 


0.0710 


-0.1415 


0.5145 


0.7448 


-0.1608 


0.0751 


20 


0.0680 


-0.1341 


0.4730 


0.7784 


-0.1589 


0.0732 


21 


0.0644 


-0.1257 


0.4311 


0.8101 


-0.1555 


0.0707 


22 


0.0603 


-0.1164 


0.3891 


0.8398 


-0.1503 


0.0674 


23 


0.0556 


-0.1063 


0.3471 


0.8674 


-0.1434 


0.0634 


24 


0.0506 


-0.0957 


0.3052 


0.8926 


-0.1348 


0.0587 


25 


0.0451 


-0.0844 


0.2638 


0.9155 


-0.1243 


0.0533 


26 


0.0393 


-0.0728 


0.2230 


0.9359 


-0.1120 


0.0473 


27 


0.0331 


-0.0608 


0.1829 


0.9536 


-0.0979 


0.0407 


28 


0.0268 


-0.0487 


0.1438 


0.9686 


-0.0819 


0.0335 


29 


0.0202 


-0.0364 


0.1058 


0.9808 


-0.0641 


0.0258 


30 


0.0135 


-0.0241 


0.0690 


0.9901 


-0.0445 


0.0176 


31 


0.0068 


-0.0120 


0.0337 


0.9965 


•0.0231 


0.0090 



55 Claims 



1. A sampled amplitude read channel for reading digital data from a sequence of discrete time sample values gener- 
ated by sampling an analog read signal from a read head positioned over a magnetic medium, comprising: 
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(a) a sampling device for sampling the analog read signal to generate the discrete time sample values; 

(b) an adaptive equalizer, responsive to the discrete time sample values, for generating equalized sample val- 
ues according to a target response; 

(c) interpolated timing recovery for generating interpolated sample values; and 

(d) a discrete time sequence detector for detecting the digital data from the interpolated sample values. 

2. The sampled amplitude read channel as recited in claim 1 . wherein: 

(a) the adaptive equalizer comprises an FIR filter comprising a plurality of filter coefficients; 

(b) the adaptive equalizer operates according to a least mean square algorithm responsive to an error value 
computed as a function of an output of the FIR filter and an estimated ideal value; 

(c) the adaptive equalizer comprises an interpolation circuit, responsive to the interpolated sample values, for 
generating the error value e^ synchronous with the discrete time sample values; and 

(d) the least mean square algorithm uses the en-or value ei^ to adjust the fitter coefficients of the FIR filter. 

3. The sampled amplitude read channel as recited in claim 2, wherein the error value ei^ is computed relative to a dif- 
ference between an output of the interpolated timing recovery and a sample value estimator responsive to the inter- 
polated sample values. 

4. The sampled amplitude read channel as recited in claim 2, wherein the error value e^^ is computed relative to a dif- 
ference between an output of the FIR filter and an output of the interpolation circuit. 

5. The sampled amplitude read channel as recited in claim 2, wherein: 

(a) the least mean square algorithm is: 

WK.i=WK-^^-XK-e,; 

(b) Wk are the filter coefficients of the FIR filter; 

(c) ^ is a predetermined gain; and 

(d) X(c is a discrete time sample value. 

6. The sampled amplitude read channel as recited in claim 2, wherein the adaptive equalizer further comprises a con- 
straint circuit for constraining a frequency response of the adaptive equalizer at a predetermined constraint fre- 
quency. 

7. The sampled amplitude read channel as recited in claim 6, wherein the constraint circuit comprises an orthogonal 
projection operation for constraining the filter coefficients relative to the predetermined constraint frequency. 

8. The sampled amplitude read channel as recited in claim 7, wherein: 

(a) the constraint circuit operates according to a modified least mean square algorithm, 

Wk.i = Wk-*^-PviV2^-{X,-eK) 

(b) are the filter coefficients; 

(c) ^ is a predetermined gain; 

(d) PviV2-*- is the orthogonal projection matrix; and 

(e) X)( is a discrete time sample value. 

9. The sampled amplitude read channel as recited in claim 1, wherein the interpolated timing recovery further gener- 
ates a data clock for clocking the discrete time detector. 

10. The sampled amplitude read channel as recited in claim 1 , wherein: 

(a) the analog read signal comprises pulses modulated by the digital data at a predetermined baud rate; and 

(b) the sampling device samples the analog read signal at a rate above the baud rate. 

1 1 . The sampled amplitude read channel as recited in claim 1 , wherein the interpolated timing recovery comprises an 
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interpolator responsive to an interpolation interval x proportional to a time difference between a channel sample 
value and a desired interpolated sample value. 

12. The sampled amplitude read channel as recited in claim 1 , wherein the analog read signal comprises pulses mod- 
ulated by the digital data at a predetermined baud rate, the interpolated timing recovery further comprising: 

(a) a frequency offset generator for generating a frequency offset Af proportional to a frequency difference 
between a sampling clock and the baud rate; and 

(b) a mod-Ts accumulator for accumulating. modulo-T$, the frequency offset Af to generate an interpolation 
interval x where Is is a predetermined sample period of the sampling clock. 

13. The sampled amplitude read channel as recited in claim 12, wherein the frequency offset generator comprises: 

(a) a phase error detector for detecting a phase en^or AO between an interpolated sample value and an esti- 
mated sample value; and 

(b) a loop filter for filtering the phase error A6 to generate the frequency offset Af, 

14. The sampled amplitude read channel as recited in claim 11, wherein the interpolator comprises an actual interpo- 
lation filter having an actual impulse response h(k} which approximates an ideal impulse response of an Ideal inter- 
polation filter: 

sinc(7c • (k-x/Ts)); where: 

k is a time index; and 

Ts is a sample period of the sampling clock. 

15. The sampled amplitude read channel as recited in claim 11, wherein the interpolator comprises an actual interpo- 
lation filter having an actual impulse response h(k) generated by minimizing a mean squared error between a fre- 
quency response of the actual interpolation filter and an ideal frequency response corresponding to an ideal 
interpolation filter. 

16. The sampled amplitude read channel as recited in claim i 1 , wherein: 

(a) the interpolator comprises an actual interpolation filter; and 

(b) a plurality of coefficients C^y of the actual interpolation filter are computed in real time as a function of x. 

17. A method for reading digital data from a sequence of discrete time sample values generated by sampling and ana- 
log read signal from a read head positioned over a magnetic medium, comprising: 

(a) sampling the analog read signal to generate the discrete time sample values; 

(b) adaptlvely equalizing the discrete time sample values to generate equalized sample values according to a 
target response; 

(c) interpolating the equalized sample values to generate interpolated sample values; and 

(d) detecting the digital data from the interpolated sample values. 

18. The method for reading digital data as recited in claim 17, wherein the step of interpolating sample values com- 
prises the step of generating a data dock for clocking the discrete time sequence detector. 



23 



EP0 805 447 A2 



USER 
DATA 



RLL 
ENCODER 



r 



PRECODER 



DATA 
GENERATOR 




SYMBOLS 
>+U a(n) 



18 

[ MAGNETIC 
RECORDING 
CHANNEL 



22 



WRITE 




58 



19 
21 



60 




GAIN 
ICONTROLJ 



50 



— M.CIRCUITRY 
20 

f ANALOG ^ ^} 
RECEIVE"-^ 

FILTER 



28 



26 

25 / 

f DISCRETE 

EQUALIZER 

FILTER 



23 



TIMING 
RECOVERY 




52 




FREQUENCY 
SYNTHESIZER 



B143 



CDR'^30 



34 ^ ^ 

W DISCRETE ^ , 
TIME I b(n) 



SEQUENCE 
DETECTOR 



RLL 
DECODER 




DATA SYNC 
DETECTOR 



ESTIMATED 
DECODED 
USER 
DATA 



FIG. 1 



PRIOR ART 



24 



EPO 805447 A2 




FIG, 2A 

PRIOR ART 



15 




FIG, 2B 

PRIOR ART 



25 



EPO 605447 A2 



2 

USER^ 
DATA* 




~b(n)<rri 



12 




SYMBOLS 
>+^^ a(n) 



MAGNETIC 
RECORDING 
CHANNEL 

22 




^21 



WRITE 
CIRCUITRY 



ANALOG , , 
RECEIVE 
FILTER 




W FREQUENCY 



50 



54 



SYNTHESIZER 



CDR'^30 




DISCRETE 
EQUALIZER 
FILTER 



r INTERPOLATED 



B102 



TIMING 
RECOVERY 



34 I 

W DISCRETE 1 , 
TIME I b(n) 



8104 



SEQUENCE 
DETECTOR 





37 



DATA SYNC 
DETECTOR 




36 



00 

ESTIMATED 
DECODED 
USER 
DATA 



F/G. J 



26 



EP0 805 447 A2 




FIG, 4A 

PRIOR ART 



27 



EP0805447 A2 



Y 



(t) 



62 



B100 



52 

f FREQUENCY 
.SYNTHESIZER 




B103 

25 , ^ 

{ DISCRETE 

EQUALIZER 
FILTER 



INTERPOLATED TIMING 
RECOVERY 

. INTERPOLATED 



DATA 



SAMPLE 



VALUES? 

25 




B167^ 



Af 



PID LOOP 
FILTER 



8161 



Y 



PHASE 
ERROR 
DETECTOR 



k+T 



ACQ/TRK 



I 



B151 




B153 



B143 



f EXPECTED 
SAMPLE 
GENERATOR 



Y 



k+T 




B141 



SLICER 



DATA :B 1 
CLOCKiP 



FIG, 4B 



28 



EP0805447 A2 




29 



EP0805447 A2 




EP0805 447A2 



U(x) 
32 




Y(x)- 



I f(i.T) I Gi(n)*U(x-n) ;ff(1.T) f(2. r) f(3.T)1 

L ^ 



i=1 n=1 B262^ 



FIG. 7 



B128 



31 



EP 0 805 447 A2 




32 



EP0805447A2 




33 



EP0805 447 A2 



DIGITAL EQUALIZER MAGNITUDE RESPONSE 











• * 

• • 

> • < 

• • 

^^^^^^^^ 


















\ 


















i 1 

1 • 1 

1 i i 


















1 1 i 

> ■ 4 


















1 \ i 

1 • 1 

! 1 \ 


















: 1 i 




















1 i 

• • 




















f i 

■ ■ 

1 i 


















■ 


1 i 











0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

NORMAUZED FREQUENCY 



F/G. 9 A 



34 



EP0805 447 A2 



DIGITAL EQUALIZER PHASE RESPONSE 





■ • • ■ 1 

• • • • i 

• • • • 1 

• • • • 1 

• • • • 1 

• • • ■ I 

• • • ■ 1 

• ■ • • 1 

• • ■ • 1 

• • * ■ 4 




• a 

• • 
a • 
a a 
) a 
a a 
i a 
f • 
t • 

» a 








V • • ■ ' 

■ ■ • • 

• X. » , , « 
» » • • « 
' ■ • • < 

» X • • 1 « 
» • • • i 

; X. • • f i 
' X. • « • • 

• • t • « 




• a 

• a 
» a 
1 « 

• a 
1 • 
» a 
I * 
1 a 
1 a 








• • • ■ 

• • « 

• • a 
» » X. • • ■ 

• • • ■ • 

• X. • • • 

• • X. • • • 

• • X • • • 

• • » • 




1 a 

> a 
• • 
t • 
1 « 
t • 
» a 
» * 
I * 
1 a 








* a ax. ■ • 
' . » • X. • • 
» » • X. • • 

a a a a 

* a a X_ a a 

* a ( X • a 
► a • X. • • 

a a X^ • • 
a a ^^a a 




■ • 

» t 
1 • 

> a 
• • 
» a 
» • 

> a 








» • a 
a , JV a 

• • »x. • 

■ a a a X. ■ 
' a a • X a 
* a a a a 
' a a a X^ a 

• a » X_ 

• a • a 
a » • X • 

• a a X a 




> * 

» a • 

> a 
» * 
» a 

> a 

> a 

> a 

> • 1 
1 • ( 

> • 1 








* a a Jt. 

a a • rX^ 
a a a 1 
a • a 1 
a a a i 

• a a 1 • 

: ; : 1 >j 


t # 
ft a 

> * < 

• < 

• ( 
a 1 

• 1 
ft • 4 
ft • 1 
^ a 







a**a««aa*a 



aaaaaaacaasaaaaai 



aaaawvaaaaaa* 
a 



a 



fta»a»a»aaa*aa 
a 
a 
a 
■ 
« 
« 
a 

a 



a 

t 



aaaaa a«aa 



a a* 



• • a* a a I 
a 
a 
t 
• 

a 



a 



0.1 0.2 0.3 0.4 0.5 0.6 0.7 

NORMALIZED FREQUENCY 



0.8 



1 



F/G. 9B 



35 



EP0805447 A2 




FIG. 10 



36 



EP0805447 A2 



25 




FIR 




INTER- 




SERIAL TO 


FILTER 




POLATOR 




PARALLEL 



0100 B122 cue 



Y(2k~1)+T 



CU8 



y., CI 52, 



01 



rif(2k-1)+T 



C174 




0144.. 

^H^^ C1 

(2k-1)+T 



Y 



2k+T 



42 



0146 
01562 

+1 ALTERNATING 

iT^01582 

i RESET 
EVERY 
5TH 



0160. 



0168 



+1 ALTERNATING 
01622 

OI680 



0176 



1 



ADAPT2i.i 



DIDIDID D 
MUX 



D D D D D 



0172 



XT 



1 i i i i 

[ MUX 



,0176. 



0174, 



2 

•ADAPT 



2j 



1 



0172, 



r/G. 11 A 



37 



EP0805 447A2 



FIR 




INTER- 




SERIAL TO 


FILTER 




POLATOR 




PARALLEL 



(2k-1)+x 



CI 78 



C100 



B122 



CI 40 



(2k.i)+T 




C106 



2k-2j-1 



0^76, 



2k- 1 




(2k.1)+T 



®2k+x 




2k+'C 



GRADIENT 
AVG 



C186 



s 



GRADIENT > 



CI 46 



C186 




FIG, 11 B 



38 



EP0805447 A2 




39 



EP 0 805 447 A2 




FIR 




INTERPOLATOR 


FILTER 






^CIOO 




B122 



0148 





0106 



+1 ALTERNATING 



RESET 
EVERY 
5TH 



ALTERNATING 



ADAPTg. 



FIG. 11 D 

40 



